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We study numerically the hole pairing induced by spin-wave exchange. The contact hole-hole inter- 
action is taken into account as well. It is assumed that antiferromagnetic order is preserved at all 
scales relevant to pairing. The strongest pairing is obtained for the d-wave symmetry of the gap. 
Dependence of the value of the gap on hole concentration and temperature is presented. For the 
critical temperature we obtain Tc ~ 100 K at the hole concentration 5 ~0.2-0.3. 
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I. INTRODUCTION 

Magnetic fluctuations are believed to be a very likely 
mechanism of pairing ii|LCjuprate superconductors. There 
have been many studiesETLl of predominantly phenomeno- 
logical nature supporting this idea. In the present work 
we study spin-wave-mediated hole pairing using results 
obtained from first principles for the undoped t-J model. 

We rbas^ our study on the results of previous 
papersJld'M It was shown in Ref. |l^ that because of spin- 
wave exchange there is an effective long-range attraction 
between two holes with opposite spins, 



C/cff(r) 



A 



A < 0. 



(1) 



In this potential there is an infinite series of two-hole 
bound states. However, they have very large sizes and 
very small binding energies and thus are not directly 
responsible for high- Tc superconductivity. Very strong 
pairing in the many-hole problem due to the same po- 
tential ([I]) was demonstrated in Ref. |Tl|, where an in- 
finite set of solutions for the superconducting gap was 
found. The strongest pairing was either in the d-wave 
or g-wave sector. The pairing induced by spin-wave ex- 
change is a long-range phenomenon. However, the attrac- 
tive potential is too singular and the wave function 
is known to collapse to the origin. On the one hand, this 
"collapse" effect substantially enhances pairing. On the 
other hand, it leads to a dependence of the superconduct- 
ing gap on short-range dynamics which cannot be studi 
ied analytically. For this reason analytical calculationsEil 
can only estimate the numerical value of the gap and 
cannot distinguish between d- and g-wave pairing (which 
have the same long-range behavior but a different short- 
range one) . In the present work, we calculate the gap nu- 



merically, taking into account both spin-wave exchange 
and contact hole-hole interaction. The d-wave pairing is 
shown to be the strongest. 

Recently, d-wave pairing was studied in Ref. |l^. Al- 
though many results are similar, we believe the spin- wave 
exchange interaction which we use is more realistic than 
the atomic limit interaction employed in Ref. |l^. 

Our paper has the following structure. In Sec. II we 
present an effective Hamiltonian of the t-J model. In 
Sec. HI we calculate the BCS-type pairing of holes at 
zero temperature. Section IV presents the results of the 
calculation of the critical temperature. Finally, our con- 
clusions are given in Sec. V. 



II. EFFECTIVE HAMILTONIAN FOR 
DRESSED HOLES 

The underlying microscopic physics is described by the 
t-J model defined by the Hamiltonian 



H 



Ht + Hj 

{nm)(T 



H.c.) -t- J ^ S„ • S™, (2) 

(nm) 



where dj^^ is the creation operator of a hole with spin 
(T ((T =t)i) 3't site n on a two-dimensional square lat- 
tice. The djj^ operators act in Hilbert space with no 
double electron occupancy. The spin operator is S„ = 
^dl^^(Ta:(3dn/3- {nm) are the nearest-neighbor sites on the 
lattice. Below we set J 1 and give all energy values in 
units of J. 

At half-filling (one hole per site) the t-J model is equiv- 
alent to the Heisenberg antiferromagnet model which 
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has long-range Neel order in the ground state. Un- 
der doping the long-range antiferromagnetic order is de- 
stroyed. However, local antiferromagnetic order is pre- 
served. We assume that the magnetic correlation length 
Cmagn IS not Smaller than the typical wavelength of holes, 
^magn > ^/pf ^ ((5 <C 1 is the concentration of 

holes) . Thus we have antiferromagnetic order at all scales 
relevant to the problem. This assumption does not con- 
tradict experimental data£3 

We treat the J term of the Hamiltonian (H) using the 
linear spin- wave approximation (see Ref. |lj for a re- 
view). Define the Fourier transformations 



net nei 

where the notation n G| [n G^) means that site n is on 
the spin-up (-down) sublattice. Introducing the Bogol- 
ubov canonical transformation 

we write the Heisenberg Hamiltonian Hj as 



Hj = £^0 + Cc)q(Q:|jaq + /3^/3q), 



q 



(5) 



where is the antiferromagnetic background energy. 
The summation over q is restricted to the Brillouin zone 
of one sublattice [7q — ^(cosga: + cos qy) > 0]. The spin- 
wave dispersion and the transformation coefficients are 
given by 



a;q = 2Jl-72 Wq w ^/2|q| at |q| < 1, 



1 1 



LUq 2 



Vq = -sgn(7q). 



1 1 



U!q 2 



(6) 



The spin waves created by and have definite values 
of spin projection. Due to Eqs. (||) and (||), a^|0) has 
= -1 and /3t|0) has 5^ = +1. Here |0) is the wave 
function of the quantum Neel state. 

Single-particle properties in the t-J model are by now 
well established (see Ref. |l^). A single hole is a mag- 
netic polaron of a small radius, i.e., a "bare" hole that is 
"dressed" by virtual spin excitations. A single hole has 
a ground state with a momentum of k = (±7r/2, ±7r/2). 
The energy is almost degenerate along the line cos + 
cos ky — which is the edge of the magnetic Brillouin 
zone (see, e.g., Refs. p^-p^). The hole dispersicffL.may be 
well approximated by the analytical expressionEiJ 



^A2/4 + 4t2(l+y) 

^A2/4 + 4t2(l+y)-4i2(x + y)7j 



+ -/32(cosfcx - COS ky) ; 
Ao « 1.33, X « 0.56, y « 0.14, 



where the parameters Ao,x,y arCpSpme combinations of 
the ground state spin correlators.Eil Near the band bot- 
tom ko = (±7r/2, ±7r/2) the dispersion (|^) can be pre- 
sented in the usual quadratic form 



P2 < (ii, 



(8) 



where pi (^2) is the projection of k — kg on the direction 
orthogonal (parallel) to the face of the magnetic Brillouin 
zone (Fig. 1). From Eq. (0) for t > Ao/4 we have 



y 



t « 0.65i. 



(9) 



According to Refs. |T| and ||, (32 « O.li at t > Ao/4. 
The wave function of a single hole may be written in 
the form ■i/'kcr = ^kcr|0)- At large t the composite hole 
operator /ij^.^ has complex structure. For example, at 
t/ J = 3, very roughly, the weight of a bare hole in ■0kcr 
is about 25%; the weight of configurations "bare hole -I- 

1 magnon" is ~50%, and of configurations "bare hole -t- 

2 or more magnons" ~25%. These estimations-are based 
on the approach with a minimal string ansatztll pad fur- 
ther renormalization due to additional magnons.o Ho"Vfe 
ever, other approaches like finite cluster diagopaliaationliZl 
or numerical solution of Dyson's equationtSllSo'Ea give 
very close results. We have to stress that the dressed hole 
is a normal fermion. 

The interaction of a composite hole with spin waves is 
of the form (see, e.g., Refs. ^8|jl9|j2|j2|) 



Hh.sw = ^ .gk,q (^k+qi^kTttq + ^k+qt^ki/3q + H.C. 



k,q 



5k,q = 2/ y -^(7kt^q -I- 7k+qV'q). 



(10) 



(7) 



For arbitrary t the coupling constant / was calculated in 
Refs. The plot of / as a function of t is presented in 

Fig. 2. For large t the coupling constant is <- independent 

Let us stress that even for t > J the interaction (nO) 
between quasiholes and spin waves has the form as for 
t ^ J (i. e., as for bare hole operators) with an added 
renormalization factor (of the order of J/t for t ^ J). 
This is a remarkable property of the t-J model which is 
due to the absence of a single-loop correction to the ver- 
tex. This property was first found perhaps in Ref. |l^. In 
Refs. p"8|Jl9|j2^ it was demonstrated explicitly that vertex 
corrections with different kinematic structure are of the 
order of few percent at t/J ^ 3. There is also a weak q- 
dependence of the coupling constant /. The plot in Fig. 2 
corresponds to the long-wavelength limit q = because, 
as we will see later, the small g's are most important for 
pairing. At g tt, the factor / is (10 — 17) % bigger than 
at g = [sec the discussion between Eqs. (13) and (14) of 
Ref. |l^ . The influence of this correction on the pairing 
is negligible. 
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Note that the hole scattering between different pockets makes a large contribution to the pairing. However, in the 
two-sublattice formalism which we use, there are no spin waves with q = g = (±7r, ±7r) and such scattering takes 
place via umklapp processes with q ^ pp <C 1. One could use another description: Expand the Brillouin zone for spin 
waves and include q « g into consideration explicitly. Then, due to antifcrromagnetic order the points q = and 
q = g are equivalent, and the coupling constants in the effective Hamiltonian (nQ) are exactly equal, /q=o — /q=g- 
Certainly the kinematic structure of the vertex (|lo| ) reflects this symmetry: gk,q = 5k,q+g- 

Interaction between two holes can be caused by exchange of one spin wave. Alongside that there is a contact 
hole-hole interaction. One can say that it is due to exchange of several hard spin-wave excitations. The Hamiltonian 
of the contact hole-hole interaction was derived in Refs. PqJlCJ using a variational approach: 



TV ^ 

ki,k2,k3,k4 



^7ki-k3 + ^C'(7ki+k3 + 7k2+k4) 



A = 16ii//i^(l - 7/1^) - - - 2/i 



2 18.5^"* + 8 V 



T ^til^^2 I ^ki T ^ki +k2 ,k3 + k4 , 



(11) 



C = -ati^fT 



where 



1 

"=2 



3/2 + 2St 
S't 



1/2 



[St{3/2 + 2St)V/^ 



St = [9/16 -I- 4i^] 



211/2 



(12) 



The coefficients A and C in Eq. ( |11| ) were derived in first order in a, where a is the coefficient in front of the transverse 
part to the Heisenberg interaction: S„Sm S^S^+^{S^ S^+S^ S^). Since the physical value is a = 1, contributions 
of higher orders are important. In order to estimate them, we will set a = 0.6. This choice is made sjj that results for 
the binding energy oLshart-range two-hole bound states obtained by finite lattice diagonalization£3"Ej would agree 
with results obtainedn3'E2l by using the effective interaction (pT|). Actually at a = 0.6 the contact interaction Hh.h is 
very small and practically does not influence the pairing. 

To summarize, we conclude that the dynamics of holes on the antifcrromagnetic background is described by the 
effective Hamiltonian 



kcr 



(13) 



which is expressed in terms of the composite hole h^a and spin-wave Qg, /3a operators. It includes free holes and spin 
waves and their interactions i/h.sw and i?h,h [given by Eqs. (^0|) and (|ll|)]. 



III. SUPERCONDUCTING STATE 



For the small concentrations S under consideration, holes are localized in momentum space in the vicinity of the 
minima of the band, ko = (±7r/2, ±7r/2), and the Fermi surface consists of ellipses (see Fig. 1). The Fermi energy 
and Fermi momentum of noninteracting holes are 

= ^7:{f3,P2)'^'6, PF ^ (7rJ)i/2. (14) 

The Fermi momentum pf is measured from the center of the corresponding ellipse. Let us stress that the numerical 
value of ep is very small. For realistic superconductors t/J fa 3 (see, e.g., Refs. ^2|"^4|). Therefore at i5 = 0.1 
and J = 0.15 eV one gets ep ~ 15 meV « 175 K. In pairing, the exchange of spin waves with typical momentum 
q ^ Pp <C 1 is the most important. The energy of such spin waves is much higher than the typical energy of a pair, 

LUg^pp^ (7r(5)i/2 » ep ~ (/3i/?2)'/'<5. (15) 

The situation is quite similar to that for the two-hole bound state problemS and much different from the situation 
with the usual phonon- induced pairing where Debye's frequency is much lower thaa.the Fermi energy. 
The interaction between two holes with opposite spins and opposite momenta isliJ 

M.,k' - -2 g^-'^y--'^ + ^(A7k-k' + C7k+k'). (16) 
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The first term here is due to the spin-wave exchange 
diagrams shown in Fig. 3. The minus sign before this 
term takes into account the fact that spin- wave exchange 
makes the spin flip for both holes. For the same reason, 
the momentum transfer is the sum (not the difference) of 
the hole momenta q = k + k'. The energy denominator 
in Eq. ( [l^ ) takes into account the energy of the spin- 
wave LOq, and the energies i?k and E'k' of the two holes 
in intermediate unpaired state. In fact, the account of 
and £^k' is the account of retardation. We discuss 
this question below. The second term in Eq. ( p^ is the 
contact interaction (|ll]). 

We use the usual BCS wave function for the ground 
state of the many-hole system 



(17) 



Thus we suppose that all quasiparticles are in the con- 
densate. For strong interactions the validity of this as- 
sumption is under question because there is no parameter 
to justify it. We believe that numerically the wave func- 
tion (|l^) is good. Anyway one may consider the wave 
function (|l7| ) as a trial one in the variational method. 
In this case the large gain in energy which we get is a 
justification of the wave function. 

The gap Ak corresponding to the wave function ([11 
satisfies the conventional BCS equation 



ye 



(18) 



where = Ck — A* being the the chemical potential 
fixed by the hole density 



(19) 



It is well known that the excitation energy of fermions in 
BCS theory is Ei^ = VCk + ^k- J^^t this energy enters 
Eq. (flq) for the effective hole-hole interaction. Equa- 
tion dig ) is obtained by variation of the average value of 
the Hamiltonian with respect to the parameters Uk and 



denominator in the spin- wave contribution is — ^^cxcited, 
and it does not depend explicitly on Wk and Wk. The 
self-consistency condition £ — fexciteci^ —^q — -Ek — -Ek' 
appears after solving Eqs. (^) and (|2l|). From the prac- 
tical point of view this question is not important because 
due to the condition ( |l5| ) the dependence of Vk,k' on the 
gap is very weak. 

An iterative numerical solution of Eq. (|l^) is straight- 
forward. We present results for Jl corresponding to 
realistic superconducting systems-Ej lj Since the inverse 
mass /?2 [see Eqs. (0),(||)] is known with rather poor ac- 
curacy, we use several values of the mass ratio a — /3i//?2- 
We take Pi from Eq. (^ and then set (32 — Pi/ a. The 
constant of the hole-magnon interaction ( p^ is / = 1.80 
at i = 3. 

The symmetry group of the square lattice is The 
solutions of Eq. (18) belong to certain representations of 
this group. In agreement with Ref. |ll|, the strongest pair- 
ing is in the Bi representation [d wave. Fig. 4(a)] and in 
the A2 representation [g wave, Fig. 4(b)]. Consider first 
the d-wave pairing. The map of the gap for the hole con- 
centration 5 — Q.\ and the mass ratio a = P1/P2 = 7 
is presented in Fig. 5. In Fig. 6 for the same parame- 
ters we give the map of which is the mean occupation 
number of a single-hole quantum state. We observe that 
despite a big value of gap the hole density distribution 
changes quite sharply at crossing the Fermi surface. For 
other mass ratios and hole concentrations Fig. 5 is also 
approximately valid because it gives the gap in units of 
Aniax and we found that with changing S and a, the whole 
gap function is multiplied by some factor but the k de- 
pendence is not much changed, 



Ak((5, a) 



c((5,a) 



7vAk(<J',a')- 



(22) 



Due to interaction between holes, the ideal gas relation 
( |l^ ) between the chemical potential /i and the hole con- 
centration is not valid. The correct relation follows from 
Eq. (p^. Plots of /X as a function of 5 are given in Fig. 7. 

A very important characteristic is the maximam value 
of the gap on the Fermi surface, Ai. Its value is directly 
related to the critical temperature of the superconducting 
transition,llil 



5u] 



-{^>\H-£\-^) 



(5uk 



0. 



(20) 



Here £ is the energy of the ground state. The effective 
interaction (|l6|) itself depends on the parameters Uk and 
Wk via the dependence of f^k on the gap Ak. Nevertheless, 
in the variational equations ( po|) we have to set 



(21) 



S S 
^Vkk' = ^Vkk' = 0, 



and therefore we get the usual BCS equation (^8|) . Expla- 
nation of the condition (|l|) is as follows. The spin- wave 
exchange part of the interaction ( p^ ) is due to the sec- 
ond order of perturbation theory. Therefore, the actual 



Tc w 0.5Ai(T = 0). 



(23) 



According to Fig. 5, Ai « 0.7Aniax- The dependence of 
Al on the concentration is given in Fig. 8. Comparing 
the plots of the gap Ai and the chemical potential fi 
(Fig. 7), we conclude that Ai ^ 0.7//. This is really 
a very strong coupling limit and virtually all holes are 
involved in pairing. This is to be contrasted with the 
usual situation when only a small portion of electrons 
^ take part in pairing and the gap is proportional to 
the Debye frequency. 

The (jf-wave pairing is weaker and we will not present 
complete results for this case. Due to the above men- 
tioned similarity of the long-range (small q) behavior of 
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the d and g waves which arises from having the same 
number of zeros at the Fermi surface, the value of the gap 
for the g wave for the above parameters is of the same 
order as for the d wave. It is also interesting that the g 
wave does not depend on details of the contact part of 
the interaction (Jl^) while the d-wave gap is substantially 
suppressed by adding repulsion to the short-range inter- 
action. Thus, under certain conditions the g-wave solu- 
tion may be relevant to the problem. Table I gives more 
information about solutions at several parameters includ- 
ing the difference of the free energy F = — /iiV;i|5') 
{Nh is the number of holes) between the superconducting 
and normal states, 

Fs~Fn = 2Y^ ^kWk + \Y1 ^k,k'Uk«kUk'Wk' • (24) 

k k,k' 

It is convenient to calculate free energy per hole and use 
the difference fs - Jn = ^{Fs - Fn). 



The most important remaining problem is the destruc- 
tion of long-range aatiferromagnetic order. Following 
experimental resultstS we have assumed that antiferro- 
magnetic order is preserved at distances r ^ i/pp- The 
behavior at larger distances is an open question in the 
present paper. 
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IV. CRITICAL TEMPERATURE 

Due to the condition fHq) the spin-wave frequency in the 
hole-hole interaction ( Jig ) is large in comparison with the 
hole excitation energy: Wq » i^k- It means that retarda- 
tion is small and the interaction is almost instantaneous. 
It is well known that in this case the equation for the gap 
at T ^ is 




In Fig. 9 we present the calculated dependence of Ai on 
temperature at hole concentration 5 = 0.1. Figure 10 
gives the dependence of the critical temperature Tc on 
hole concentration. The approximate relation ( |23|) de- 
rived analytically in Ref. |ll| is qualitatively fulfilled. In 
real units (J = 0.15 eV), Fig. 10 gives (taking a — 7) 
Tc = 51 K at (5 = 0.1 and = 86 K at ^ = 0.3. Let us 
stress that in our calculation we do not use any fit. The 
only input is the values of t and J. 

V. CONCLUSIONS 

Using the single spin-wave exchange mechanism sug- 
gested in Refs. ^|,^we carried out a numerical ab initio 
calculation of superconducting pairing in the t-J model. 
Both the magnitude of critical temperature and its de- 
pendence on hole concentration are in good agreement 
with experimental data. The calculated critical temper- 
ature is still smaller than the highest critical temperature 
obtained in experiment. However, this may be explained 
by not knowing the exact parameters. By a relatively 
small variation of parameters we can get Tc = 100-150 
K. 
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TABLE L Influence of short-range interaction. Changes in 
the short-range interaction are introduced by increasing the 
parameter a [see discussion after Eq. (^l])] which makes the 
contact interaction more repulsive. All data are for T = 0, 
a = 7, fi = 0.058 (this is the chemical potential at which in 
the absence of interaction the hole concentration would be 
So = 0.05). 



FIG. 1. The Brillouin zone of a hole in the t-J model. 
FIG. 2. The plot of the hole-spin- wave coupling con- 
stant /. 

FIG. 3. Interaction between two holes via a single spin- 
wave exchange. 

FIG. 4. The gap symmetry, (a) Bi type {d wave), (b) 
A2 type {g wave). 

FIG. 5. The contour plot of the d-wave gap for t/ J ~ 3, 
the mass ratio a = /3i//32 = 7, and hole concentration 6 — 
0.1. The levels are presented in units of the gap maximal 
value which at these parameters is equal A^ax = 0.0661. 
Dashed curves represent the Fermi surface for the case 
when one considers the holes like an ideal gas. 

FIG. 6. The contour plot of a single hole quantum 
state mean occupation number v^. Parameters are the 
same as in Fig. 5. 

FIG. 7. The chemical potential /i as a function of hole 
concentration 5. Dashed curves correspond to an ideal 
gas of holes. Deviation of dashed curves from linear de- 
pendence (|l^) is due to the deviation of the dispersion 
relation (^) from the quadratic expansion (H). The de- 
pendence with pairing taken into account is presented by 
solid lines. All the curves correspond to t/J = 3. The 
mass ratio is (from top to bottom) a = P1/P2 = 5, 7, 9. 

FIG. 8. The maximal value of the gap on the Fermi 
surface Ai vs hole concentration for t/J = 3; the mass 
ratio is (from top to bottom) a = /3i//32 = 5, 7, 9. 

FIG. 9. The maximal value of the gap on the Fermi 
surface Ai as a function of temperature. The hole 
concentration is (5 = 0.05, t/J = 3, the mass ratio 
a = /3i//32 = 7. 

FIG. 10. The critical temperature vs hole concentra- 
tion. t/J — 3, the mass ratio a = Pi/ 02 — 7. 
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